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A decrease in the rotational period observed in torsional oscillator measurements was recently taken as a 
possible indication of a putative supersolid state of helium. We reexamine this interpretation and note that the 
decrease in the rotation period is also consistent with a solidification of a small liquidlike component into a low- 
temperature glass. Such a solidification may occur by a low-temperature quench of topological defects (e.g., 
grain boundaries or dislocations) which we examined in an earlier work. The low-temperature glass can account 
for not only a monotonic decrease in the rotation period as the temperature is lowered but also explains the peak 
in the dissipation occurring near the transition point. Unlike the non-classical rotational inertia scenario, which 
depends on the supersolid fraction, the dependence of the rotational period on external parameters, e.g., the 
oscillator velocity, provides an alternate interpretation of the oscillator experiments. 

PACS numbers: 73.21.-b 



I. INTRODUCTION 

Supersolidity is a unique state of matter that simultaneously 
displays both superfluidity and crystalline order. Much of our 
understanding of this exotic state is based on the pioneer- 
ing work of Penrose and Onsager^i Andreev and Lifshitz^ 
Chester and ReattOiMaS Leggett^^ and Andersonj^ Recently, 
Anderson and coworkers revisited this problemi^ Current de- 
velopments in this field, presented at a KITP workshop at 
Santa Barbara, are available onlinCiiS 

*He has long been thought to be the most likely candi- 
date for the supersolid state. Torsional oscillator experi- 
ments by Kim and Cha n"''~'^''^ generated renewed inter- 
est in this possibilityji^ In addition to the work of Kim 
and Chan j"''"''^-'"* there is now an independent confirma- 
tion of the anomalous behavior of solid ^He, as presented by 
the groups of Reppy and Shirahama i ' ' 1^.20.2 1 ^jj Qf these 
groups use torsional oscillators similar to those employed by 
Kim and Chan i'^-'^'^^-'"^ Rittner and Repp}^ report history de- 
pendence of the signal, when annealing the sample, to the 
extent of no observation of any mass decoupling in the tor- 
sional oscillator experiment. On the other hand, rapid freez- 
ing of helium leads to disorder and a drastic increase of the 
signaL™ Taken at face value, these torsional oscillator exper- 
iments indicate an anomalous mechanical behavior of solid 
^He at low temperatures. Nevertheless, the connection be- 
tween the mechanical measurements and the suggested super- 
solidity remains tenuous. The case for supersolidity entails a 
certain assumption, which we will discuss at length below, for 
the interpretation of the oscillator data. 

The measurements by Refs. [Tllflllll4lflMll20ll2ll mon- 
itor the period of torsional oscillators in the presence of an 
applied torque. The purpose of our work is to examine the 
mechanical properties of the torsional oscillator. We will only 
allude to linear response theory and causality (the Kramers- 
Kronig relations). We show that the change in torsional os- 
cillator period may be triggered by a change in damping and 
other oscillator parameters. Such a change may result from a 
low-temperature quench of residual topological defects in the 
solidification into a glass. Here, we will use the term "liquid- 



like" component to describe elastic defects that are not pinned 
above a glass transition at high temperatures. In Ref. |22[ 
we first explained how a low-temperature pinning of elastic 
defects (two-level tunneling systems) can account for the re- 
ported excess low-temperature specific heat with regards to a 
perfect Debye crystal. The results of Ref. ^ are at odds with 
a simple uniform supersolid transition, but can also occur in 
an exotic glassy supersolid state as discussed by Boninsegni 
et al^ or Wu and Phillips^l. The effects described in this 
paper may account for the observed change of the torsional 
oscillator period. The central thesis of our current work is that 
there is no imperative to adduce a supersolid fraction which 
depends very sensitively on a host of external parameters in 
order to explain the reduction in the torsional oscillator pe- 
riod as temperature is lowered. Rather, the observed trends 
might be a very natural and quite universal outcome of a glass 
transition (regardless of its classical or quantum nature). 

A direct proof of superfluidity would be the observation of 
persistent current. In this regard, we note a recent experimen- 
tal search for superflow by Day et. al.r^j^ who found no mass 
flow of any kind to very high accuracy. Relying on these re- 
sults and, most notably, on a new thermodynamic analysis, 
we recently concluded?- that the effect first observe d"-'^-'^ 
is likely not an intrinsic property of solid ''He. This is so 
because the low-temperature behavior depends critically on 
the ■^He concentration, shows annealing and rapid freezing 
dependence^iSi^i and, most importantly, has a specific heat 
which is inconsistent with a simple supersolid transition of 
a 1 % condensate fractioni22 However, our earlier work2^ does 
not rule out a supersolid fraction of order 0(10"^). 



II. OUTLINE 

Our principal findings are briefly summarized in Sec. dlllb . 
To set the stage for these results, we will briefly review, in 
Sec. dlVb . the non-classical rotation inertia (NCR!) effect. We 
then discuss the cuiTent possib le grou nd-brea king irn plication 
to the observations of Refs. [TTIIiai3lfl4tfl8lfl9| 2oi2l| [wherein 
the change in torsional oscillator period may be triggered by 
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a supersolid phase]. We then step back in Sec. (|V]l and an- 
alyze the resuhs for a first principle linear response analysis 
without assuming the NCRI. Consequently, we show how the 
experimental results are consistent with a change in damping 
or other oscillator parameters which are brought about by the 
solidification of a small liquidlike component as the temper- 
ature is lowered. In Sec. dVII i. we examine the experimen- 
tal consequences of a uniform supersolid transition and illus- 
trate that a simple supersolid transition cannot account for the 
peak in the dissipation observed by Ref. [l^ We then deter- 
mine, in Sec. dVIII ). the angular response for a glass transi- 
tion, which can account for all of the existing trends in the 
data, derive in detail our new results on the dissipation and 
torsional oscillator period for a glass, and suggest an experi- 
ment to distinguish between the glass versus supersolid sce- 
nario. We conclude in Sec. (IVIIIb with a brief summary of 
our results and their relevance to the torsional oscillator data 
of Refs. [TllfTlTll4lfl8lfl9ll20l^ In the appendix (Sec. \B, 
we illustrate that other, non-glass, response functions can also 
account for the data without assuming a supersolid transition. 



in. PRINCIPAL RESULTS 

We will address the linear response theory of a driven tor- 
sional oscillator and derive our key results in Sec. ( I VIII ). A 
general treatment of the problem involves the analysis of pos- 
sible modes in the vessel which are excited by the external 
torque and their influence on the torsional oscillator. Depend- 
ing on their character, such modes can easily be disrupted by 
the insertion of barriers and lead to a far smaller difference 
between the response of a liquidlike component and that of a 
glasslike component. To make the discussion clear, we will 
focus on a simple form for these modes, which are character- 
istic of glass formers. Here we will show: 

A glass transition, irrespective of whether or not it appears 
in a normal solid or supersolid, leads to specific testable pre- 
dictions for the oscillator period and dissipation.-^^ -'* Unlike 
a simple supersolid transition, i.e., a transition where only a 
change in the inertia of the oscillator occurs, our key result 
Eq. (|24] | predicts that, as the temperature is lowered, two prin- 
cipal trends will be noted: (i) a monotonic decrease in the 
oscillator period and (ii) a peak in the dissipation. A uniform, 
mean-field type, supersolid transition cannot account for phe- 
nomenon (ii). 

It is important to stress, that unlike all of the works to date, 
we do not assume that only a supersolid can account for the 
observed decrease in the oscillator period [feature (i)]. In- 
deed, as we show for the first time, glassy relaxations man- 
date the current observations of features (i) and (ii). This 
occurs universally and is independent of the specific micro- 
scopic mechanism by which relaxation processes occur or the 
detailed character of the glass. It is nevertheless encouraging 
to note that samples in which the most pronounced anomalous 
behavior is observed contain defectsi^ii^i 

We considered other (non-glass) "freezing" transitions that 
can also account for the empirical trends. To make this con- 
crete, in the appendix A, we will detail how a decrease in 



damping (as a result of a freezing transition) can also account 
for the observed tendencies. Unlike for the results discussed 
above, here we cannot provide a specific functional form to fit 
to. To provide the simplest qualitative form, we invoke Gaus- 
sian distributions for the quantities in question. 

Finally, we note that throughout our work, our general form 
for the angular response allows for a dependence on external 
parameters such as the initial oscillator velocity. These param- 
eters, e.g. stiffness and density and cell torsional frequency, 
which parameterize the external torques, are convolved with 
the intrinsic angular response. A liquidlike component in the 
cell may respond differently to different initial oscillator ve- 
locities. Such a dependence is observed experimentally. Aside 
from other possible effects, external parameters [whether in a 
supersolid, a glass, or in any other phase] will also always 
appear via their incorporation in the external torques in the 
general form for the angular response. 

IV. NON-CLASSICAL ROTATIONAL INERTIA IN A 
NUTSHELL 

The basic idea underlying the torsional oscillator measure- 
ments is a test of the non-classical rotational inertia (NCRI) 
- the inability of a superfluid to rotate under an applied tor- 
sional drive. The analog of this experiment on liquid ^He was 
performed by Andronikashviliis^ Later Hess and Fairbanks 
verified the existence of the NCRI and of a rotational Meiss- 
ner effect for superfluid helium, which was predicted by 
London.— As a consequence of a superfluid state, which no 
longer rotates with the "normal" liquid, the effective moment 
of inertia decreases as the superfluid fraction fg increases 

leff ^ Iosc+lHe[l- fs{T)] (1) 

which is smaller than the "classical" moment of intertia / = 
lose + I He of the combined oscillator-helium system. For an 
ideal torsional oscillator of stiffness a, a measure of the period 
(P = 2T:yJleff /a) indicates what the superfluid fraction fs 
is. The inferred results for fj^ agree well with the superfluid 
fraction measured by other probes. The current situation for 
solid "'He is far from being as clear cut. The experiments show 
that the rotational motion under an applied torque becomes 
more rapid at lower temperature- suggesting an analogue of 
the NCRI for a supersolid phase of helium. In this paper we 
point out that the measurements are also consistent with less 
exotic material effects, which can be distinguished by further 
experiments. 

V. HOW CAN WE SIMPLY UNDERSTAND THE CHANGE 
IN TORSIONAL OSCILLATOR PERIOD IF IT IS NOT 
ASSOCIATED WITH A SUPERSOLID? 

The NCRI interpretation of the existing data is very sim- 
ple and suggestive. In what follows, we present a very brief 
mechanical analysis of the measurements which points to an 
equally simple explanation. To keep the analysis as general 
as possible, we rely on first principles alone. The torsional 
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oscillator experiments measure the susceptibility - they do not 
directly monitor the moment of inertia of the supersolid. We 
start by writing down the general equations of motion for a 
torsional oscillator defined by an angular coordinate 9, 



\losc ~ 



a 



W) = 



Tecctit)+ [ ~g{t,t'-T)e{t')dt'. 



(2) 



Here, lose is the moment of inertia of the torsional oscillator, 
a is its restoring constant, ^/osc is the dissipative coefficient of 
the oscillator, g arises from the back action of solid ''He (it 
plays the role of a "polarization"), and Text is the externally 
imposed torque. In general, g is temperature (T) dependent. 
The oscillator angular coordinate 9{t) is a convolution of the 
applied external torque ff.x and the response function x{t, t'), 



~6{t) ^ / dt' x{t,t')Text{t'). 



(3) 



Causality demands that x{t < t') ~ 0. This implies that 
x{t, t') = 9{t—t')x{t, t'); under a Fourier transformation, this 
leads to the Kramers-Kronig relations which we will briefly 
touch on later. As in any time translationally invariant system, 
the Fourier amplitude of the angular response of the torsion 
oscillator is a product 



(uj) = xi^^)Textiuj), 



(4) 



with X ~ Xi + *X2 the angular susceptibility and Text the 
external torque in Fourier space. For the simple torsional os- 
cillator. 



X 



(5) 



Here, g{u!,T) is the Fourier transform of the real time g of 
Eq. (|2]). For an ideal normal solid, with a moment of inertia 
equal to Ins, which rotates with the oscillator, the back action 

ff(^;T)is 



gss{l^]T) = i"tHel^ + Ins{T)uj'^ 



(6) 



with Ins {T) the normal component of solid ^He, which varies 
with temperature and the dissipation ^He is constant. The mo- 
ment of inertia contribution of the helium is orders of magni- 
tude smaller than that of the empty vessel - the changes ob- 
served in the period are very small and are a remarkable ex- 
perimental triumph. The subscript in gss refers to the super- 
solid interpretation of the results - only a fraction of solid he- 
lium (the "normal part") rotates in unison with the oscillator 
Eq. (|6]l is consistent with the NCRl form of Eq. ([T]i- However, 
we do not need to impose this form on the existing data. Cur- 
rent experiments measure the oscillator period P = 2tt/ujo 
with ujQ the real part of the solution of 



(7) 



at fixed Ti22 For example, a decrease in the dissipative com- 
ponent jHe in the back action g{oj; T) with 



as the temperature is lowered, will also lead to a shorter rota- 
tion period. 

This observation offers a qualitatively different interpre- 
tation of the experimental result s "•'^''^•'^''^''^•^°-^' with no 
need to invoke supersolidity. Note that no superfluid phase 
is needed to explain the data and there is no NCRI effect; the 
moment of inertia Ins is temperature independent. 

The origin of the subscript in Eq. ([8]) alludes to a scenario 
wherein a higher temperature mobile liquidlike component of 
the sample "freezes" at lower temperatures leading to a de- 
crease in the dissipation 7. In the appendix, we examine in 
detail the consequences of Eq. (O. To make our discussion 
more lucid, we remark that a qualitatively similar effect ap- 
pears in the rotation of hard versus soft boiled eggs. The solid 
hard boiled egg rotates faster than the, liquidlike, soft boiled 
egg- 

The measured decrease in the rotation period 
foun d 1 1 1 1 2, 1 3. 14, 1 8, 1 9.20.21 Qjjjy implies a crossover in x 

(and a constraint on g). As the real and imaginary parts 
of X are related by the well known Kramers-Kronig (KK) 
relations. 



2 /'°° 
Xi(^) -P i du 
Jo 

X2{UJ) = ^—P I div' 



I ^'X2(w') 



a;'2 - u;2 
XI (^') 



(9) 



an enhanced decrease in xi('-^) often appears with a pro- 
nounced peak in X2(w)- A nonvanishing X2(w) at finite fre- 
quency mandates dissipation. We expect the data to indicate 
an increase in dissipation concurrent with the reduction of tor- 
sional oscillator period. Whether or not experimental data ad- 
here to the KK relations is a powerful check that needs to be 
done. 

In the sections that follow, we explain in some depth how 
the supersolid and glass pictures can both account for the de- 
crease in the rotational period when the temperature is low- 
ered. We further examine the dissipation. 



VI. A SUPERSOLID ORIGIN 

We now determine the experimental consequences of a sim- 
ple supersolid transition in which the moment of inertia fol- 
lows Eq. ([T]) (the well known NCRI effect). Here, the response 
is given by Eq. We show that while a mean-field type 
NCRI effect can, as is well known, account for the decrease 
of the rotation period it cannot account for the peak in the dis- 
sipation Q^^ observed by Ref.fisl 

We start by reviewing results for an underdamped harmonic 
oscillator. The temperature dependent period of a general 
damped oscillator of combined effective moment of inertia 



In 
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Stiffness a, and dissipation 7 is 



6{t) = i?e{Aexp[-iwo< - nt]} 



(10) 



gdissiuj; T) = i-iHe{T)uj + /„sW^ 



(8) with A a complex amplitude. For an underdamped oscillator. 
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the period 



P: 



27r 

LOq 



V4a/.//(T)-72(r) 



and the damping rate of the oscillation amplitude 



1 



21, 



(11) 



(12) 



The "quality factor" Q, which monitors the number of oscilla- 
tions required for a system to have its energy drop by a factor 
of e^^, is 



7 



(13) 



We can rewrite the changes in oscillator period as result of 
small changes in 7, /g/ / as 



SP 



An 



S-Kleffa 



+ 



4:71-^ leff 



{4.aleff - 72)3/' 



-Sa 



(14) 



We see from Eq. (fl4] i. for damping 7 < ^j2P^JJct, that the 
observed decrease in the rotation period, i.e., bP < 0] as the 
temperature is decreased may be explained by either (1) a pur- 
ported supersolid transition with g = gss, where Slejf < 
and Sj ~ 0; (2) a (non-glass) freezing transition with a sim- 
plified g ~ gdiss, where -fHe becomes smaller as the tempera- 
ture is decreased, i.e., ^7 < and 6Ief / = 0, for more details 
see the appendix|Al (3) a glass transition with an effective in- 
crease in the stiffness constant, where da > as temperature 
is lowered. The glass response to be studied later in Sec. IVIII 
emulates, in part, the features of (3). 

First, we will consider a supersolid in which g ~ gss of 
Eq. ©J 7 = 7osc + iHe is fixed and lef / (T) varies with tem- 
perature. Later on, we will comment on the effect of a drop 
in dissipation expected to accompany the supersolid state. For 
g = gss, we can deduce leff (T) and then predict the damping 
rate Kss{T). In this case, the damping rate 

lose + iHe 



2{Iosc+Ins(T)) 



aP^ + ^(aP2)2 _4^2^2p2 



(15) 



or, equivalently. 



5;2p2 + ^a4p4 _ 4,^2^,2^2 p2 



(16) 



In Eqns. ( [TSlfTSI l, only the period P varies with temperature. 
In this simple supersolid picture, the Q factor is monotonic in 
temperature. If /g^j decreases as the temperature is lowered, 
then the Q factor will decrease monotonically as the tempera- 
ture is lowered. 



A simple supersolid transition is inconsistent with the 
datai^ where a pronounced peak in Q^^ is seen near the tran- 
sition. A concurrent monotonic drop in the dissipation is ex- 
pected to accompany the supersolid transition. If the drop in 
the dissipation 7 is proportional to the drop in the effective 
moment of inertia then Q^^ will still be monotonic. Here, 
may increase as T is lowered. Such an amended simple 
supersolid picture is still inconsistent with the observed peak 
of Q^^.'** Obviously, a nonuniform change in both 7(T) and 
Ip f f (T) can be engineered to account for the data. Similarly, 
a broadening of the transition parameters, similar to that dis- 
cussed in the appendix, can lead to a peak in Q^^. 



VII. A GLASS ORIGIN 

We next illustrate how a glass transition can account for the 
observed data. A glass transition will not only give rise to a 
decrease in the period but, unlike the supersolid scenario, will 
also lead to a peak in the dissipation observed by Ref. [l^ 

The external torque is the derivative of the total angular mo- 
mentum. 



(17) 



where r is the distance to the axis of rotation, p{x) is the mass 
density and ^d{x) the local angular velocity about the axis 
of rotation. The experimentally measured quantity is the an- 
gular motion of the torsional oscillator - not that of the bulk 
helium, which is enclosed in it. Ab initio, we cannot assume 
that the medium moves as one rigid body. Similar to a ves- 
sel partially filled with a liquid component, e.g. a soft boiled 
egg, an initial imparted external torque can lead to a differen- 
tial rotation between the outer torsional oscillator and a liquid 
within it. If the liquid "freezes" into a glass the medium will 
move with greater uniformity and speed. This may lead to an 
effect similar to that of the NCRI, although its origin is com- 
pletely different. Furthermore, as a function of initial rotation 
speed, the variance between the response of "hard" glasslike 
and "soft" liquidlike media changes as temperature is lowered. 
Such a difference is indeed observed in the torsional oscilla- 
tor experiments on ^He. Our explanations do not require the 
supersolid fraction to depend on the rotation speed in order to 
account for the data. 

In what follows, we will analyze the effective equation of 
motion for the torsional oscillator. Excited modes within the 
medium in its liquid phase can lead to additional dissipative 
torques acting on the torsional oscillator. These modes can 
become far slower and effectively disappear as the system 
freezes into a glass state. Depending on the character of these 
modes, the insertion of a barrier may or may not lead to a 
change in the angular motion of the torsional oscillator. For 
modes corresponding to internal flows around the axis of ro- 
tation, such barriers will disturb the liquidlike modes. Such 
constrictions will lead to a far smaller difference between the 
period of rotation of the glass and the period of rotation with 
liquidlike component of '*He, if the modes are local. 
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Henceforth, we will focus on conventional forms for 
glasslike relaxations which might be of use for fitting data. 
Glasses are characterized by universal response functions over 
a wide range of frequencies and temperatures. We assume that 
the additional dissipative modes triggered by rotation have a 
similar form. To be specific, we will assume that g{uj;T) 
obeys the simple dependence as seen in the far better mea- 
sured dielectric response functions (where it reflects the po- 
larization) of glasses. It is given by a Dyson-like equation, 
which has been invoked in Eq. (|5]l. 



X = Xq + Xa9giXQ + ■■ 

-1 



Here, 



Xo' 



Xo 



9gi- 



[a - i-fuj - leff 



(19) 



is the inverse susceptibility of the underdamped system 
formed by the torsional oscillator chassis and normal solid, 
while ggi is the response function of the overdamped glass 
component. Physically, the first term (xo) in Eq. (fTsl l cor- 
responds to the response of the bare oscillator to an applied 
torque, the possible excitation (by the rotating oscillator) on 
a transient mode in the liquid/glass medium, which then acts 
back on the torsional oscillator at a later time (xoffgZXo); and 
all similar higher order processes. Employing the standard 
overdamped form for glass response functions^^i^^*^ 



fjgl = .90 [1 - it^s] 



(20) 



we will be able to provide an explicit form for ™d thus 
for all periods and dissipation as a function of temperature. 

In Eq. (|20] |. s is the characteristic equilibration time of the 
liquidlike component when it is perturbed. This time s di- 
verges in the glass phase (at temperatures T < Tq) and is 
nearly vanishing at very high temperatures. Similarly, < 
/3 < 1 is an exponent determining the distribution of local re- 
laxation processes. /3 corresponds to a stretched exponential 
exponent in real time. For smaller value of (3, the distribution 
of local relaxation frequencies about the dominant process at 
u!g = 1/s is more smeared out. If /3 = 1, then we will have 



a single overdamped oscillator of damping time s- Eq. 
then corresponds to the response function of an oscillator (e.g. 
Eq. ( fT9] l) with no quadratic-in-oj (inertia) term. 

The mobile liquidlike component of the sample "freezes" 
at lower temperatures into a solid glass. This, e.g., can be 
triggered by the quenching of dislocations or other defects at 
low temperatures. In Ref.22i marked superfluid-like behavior 
was indeed observed only in samples which contained grain 
boundaries - no superfluid behavior was noted in crystals of 
high quality. 

This inverse susceptibility of the glass can be approximated 
by a stretched set of overdamped oscillators. This is the physi- 
cal content of Davidson-Cole form of Eq. ( l20l i. The Davidson- 
Cole for m^^-^^'^^ fits the dielectric response data of glasses 
very well. More sophisticated nearly universal response func- 
tions, which may improve on the analysis that we give here, 
were reported in Refs. [jtIIssI In what follows, we will em- 
bark on an initial Davidson-Cole analysis of the "'He data. If 



we invoke the Davidson-Cola^i^^*^ form for ggi (Eq. ( |20] |) in 
the sum of all possible intermediate contributions of the glass 
to the total response of the oscillator, then, from Eq. ( flSl l. 



X ^i^) = [a-iju- 



.90 



(1 - ius)!^' 



(21) 



It is easy to verify, by employing Eqns. ( fTTl [T2l [T3] ), that 
for very high T (where s is small) Eq. dTTT i predicts a period 



P = 



(18) andaQ-factor 



Q 



7 + go/3s 



(22) 



(23) 



Similarly, at very low T (where s diverges), the corresponding 
quantities are given by Eqns. ( fTTl [T2l \T3[ with no change in 
the parameters given therein. A comparison of the two limits 
of high (s — ^ 0) and low (s 00) temperatures reveals that 
the Q factor is higher at low T while the period P is higher at 
very high temperatures T. 

At intermediate temperatures (where s is of the order of the 
period of the underdamped oscillator), the angle is a super- 
position of a many modes stretched about oj = 1/s. From 
Eqns. (HEI, 

1 



50 



9(t) — — / duj [a — ijuj — leff^ 

ZTT ' 



dt' Te^t{t') exp[~iuj{t - t')]. (24) 



Eq. iTM represents the final form for the angular evolution 
for all T. Here, Text{t') is the value of the external torque 
at time t' . For a delta function torque, T^xtit') = A5{t'), 
the last integral simplifies. Within our general formulation, 
and in Eq. (l24l i in particular, the observed dependence of the 
period on the external parameters is perhaps not as surpris- 
ing as demanding a superfluid fraction strongly depend on the 
oscillation velocity: For example, different initial velocities 
would correspond to different values for the internal integral 
/ dt' Text{t') cxp[—iuj{t — t')]. This integral would multiply 
the total susceptibility in all cases [whether we have a glass 
transition, a supersolid transition, or any other transitions]. In 
a treatment which is more detailed than the one which we out- 
line here, the transient modes and the consequent form of the 
response function g would also depend on external parameters 
such as initial rotation velocity. 

We anticipate that when the relaxation time is similar to the 
period of the underdamped oscillator, the dissipation will be 
maximal. This will lead to the largest magnitude in the imag- 
inary component in the singularities of x~^- Here, the glassy 
components respond with the same frequency as the "normal" 
component. At both much higher or much lower tempera- 
tures, they merely change the net effective spring constant but 
do not lead to additional modes which closely interfere with 
the oscillations of the "normal" (non-glass) part. This leads to 
a strong decoherence and to a sizable energy loss. 





FIG. 1 : (Online color) The dissipation (inverse of the quality factor) 
is shown as a function of temperature for solid ''He at 27 bar. The 
circles are data from Rittner and Reppy^ and the line is our fit 
above T > 0.1 K for a simple glass model, Eqns. | |25II26| |. with the 

functional form = As[l + (27r/os)^]^' + Q^. The fit param- ^j^-i _|_ (27r/os)'^] 



eters are A — 0.02776 sec~^, oscillator frequency fo = 184.2313 
Hz, asymptotic inverse quality factor = 11.055 ■ 10~^, glass 
relaxation time s = 1.3537e^'"^^"^ /xs, and activation barrier 
A/ks = 0.9202 K. 



We can examine the T < Tq (the glass phase) and the 
T = oo limits and perturbations about them when no dissipa- 
tion is present- 7 = 0. In these limits, we have ideal harmonic 
oscillators with no dissipation and consequently Q^^ = 0. 
When we expand about the T = 00 limit and allow for a 
small s, then as seen from Eq. ( l23b . Q^^ becomes small but 
finite (nonzero). Similarly, when T = Tq'", where s is ex- 
tremely large but finite, x has singularities off the real uj axis. 
In Eq. (|24] |. this small imaginary component now leads to a 
smaU dissipation. 



FIG. 2: (Online color) The relaxation period is shown as a func- 
tion of temperature for solid ''He at 27 bar. The circles are data 
from Rittner and Reppyi^ and the line is our simple glass model, 
Eq. ini . with the functional form of Eq. l l28t . P = {fo — 

Here the only adjustable parameter is 



B = 0.57 sec ^, otherwise same parameters as in Fig.[T] 



This Lorentzian peaks, as in Fig.[T] when s = I/luq- At this 
value of s, the poles have a maximal imaginary component 
and the dissipation is the largest. In the real dissipative oscil- 
lator (with small 7 > 0), the inverse Q value is given by 



-1 =2|xU/^ 



(26) 



with Q 



-1 _ 7. 



the linear term in 7 (see, e.g., Eq. (fT3])). 



By contrast, the correction to the absolute values of the real 



30 r 



This can be made precise for (3^1 where the poles are 
determined by the roots of a cubic equation in oj. Here, we 
can simply look for the largest magnitude of the imaginary 
part of all of the poles and see when it is maximal as a func- 
tion of s. A larger imaginary part implies a shorter decay time 
and a smaller value of Q^^. The exact expressions are not 
illuminating. The asymptotic corrections are, however, very 
transparent. As can be seen from Eq. ( |2T] ). for an ideal dis- 
sipationless oscillator (7 = 0), ojq = y/cn/ lef / is the pole 
of xi'^) in the low-temperature limit where s 00. If we 
expand about this pole, u; = ujq + ix + y, then we will find 
that the two dominant poles attain an imaginary component of 
magnitude 



90 



21, 



1 + (a;os)2 



(25) 



Rittner & Reppy 

_ theory (f =184.2313 Hz) 




10 15 



20 



FIG. 3: (Online color) The dissipation-period plot at fixed frequency 
fo — P^^^■ The curves of Figs. [T] and [2] are reploted against each 
other. The skewness of the experimental data (circles) compared to 
the semi-circle line (theory) is a well-known consequence of a real 
glass exponent /3 < 1 [see Eqs.l|20||2Ttl. which for simplicity was 
set to /3 = 1 in our analysis. 
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part of the poles is 



y 



90 



1 + {ujos)^ 



(27) 



The function y is a monotonic function of s and thus of tem- 
perature. The period 



P = 



2n 



^o + y 



(28) 



decreases monotonically when the temperature is lowered, see 

Fig.m 

The correction for the period due to dissipation is quadratic 
in 7. This will lead to (i) a monotonic decrease in the period 
as the temperature is lowered [the increase of y with lower T 
(higher s)] concurrent with (ii) a peak in the dissipation Q^^ 
at a temperature T*, for which s(T*) = I/loq (seen from the 
Lorentzian peak for |a:;|), and not to be confused with the ideal 
glass transition Tq. The inclusion of finite dissipation (7 7^ 0) 
does not modify these trends. Generally, in glasses, s follows 
the Vogel-Fulcher-Tamman (VFT) form 



s~ socxp[i?ro/(T-ro)] (r>To) 



(29) 



Here, Tq is the temperature at which an ideal glass transition 
occurs. At temperatures T < Tq, the relaxation time is as- 
sumed to be infinite. In some glasses, s may diverge alge- 
braically in (T— To). It is important to stress that, in principle, 
this may even occur if there is no finite temperature glass tran- 
sition at all and we have a simple activated form with Tq — 
in Eq. 1 



s = So cxp[A/fcBT], 



(30) 



with A the activation barrier which replaces the constant 
{DTo) of Eq. dig. In Figs. [U and |2] we plot the results of 
our perturbative calculation, Eqns. (l26ll27ll28T l. along with the 
activated form of Eq. ( [30l l and compare these to the measured 
data of Ref. 18, In Fig. [3] we show the dissipation-period plot 
of the same data at fixed frequency /o instead of the conven- 
tional Cole plot^^i^^i^ because experimentally neither xi (w), 
nor X2(w) are available at fixed T. Note that for a real glass 
there are many more parameters than those that we fit with 
here. These parameters are (i) < /3 < 1, which we set to 1 
in our equations, and (ii) the ideal glass transition temperature 
Tq > 0, which we set to zero in the activated form of Eq. dSOl l. 
In fact, an exponent /? < 1 can explain the skewness^ii^^i^ of 
the experimental result in Fig. [3] compared to the semi-circle 
result of our simplified analysis with (3^1 [see Eqs.( |25l[26l 
127] 121. 

The four relations of Eqns. (fTTl [T3l l22l l23T l along with the 
four measured values [the periods at the high and low tem- 
peratures, the Q factors at high and low temperatures] allow 
us to determine the four parameters a, 7, leff, and g. In fact, 
due to the scale invariance of the relations, there are only three 
unknowns (e.g. a/ I^f f / lef f, and g/I^ff) from which the 
fourth value can be predicted and compared to its experimen- 
tal value. Currently, the available data centers on the region of 
the transition and it is not certain what the asymptotic low and 



T+ 



high temperature values of the Q factors (Qo for T 
Qoo for T > To) or the periods (Pq for T = T+, P^ for 
T S> To) are. If we solve for all unknown quantities in terms 

of {Pq, Poo, Qoo}, then wc havc 



7//e// = MQooPoo)-'{^-Q^'r'^'' 



2ti/{QooPo. 



(31) 



The choice for the numerical value of the effective moment 
of inertia I^f / corresponds to a different choice of units. Nei- 
ther the precise values of the parameters so. A, and To of 
Eq. ( [29] l, nor of /3 in Eq. (|24] |. change qualitatively the form 
of the fits that result from Eq. (|24] |. All that matters is that 
at the temperature of the apparent transition T*, the dominant 
relaxation time 



s(r*) = 27rs(T*) ~ Fo, 



(32) 



matches with the torsional oscillator period. Putting all of the 
pieces together, Eq. (l24l i leads to a monotonic decrease of the 
rotation period as the temperature is lowered. Hand in hand 
with that, the dissipation becomes maximal when s is equal to 
the period of the unperturbed oscillator 

A possible testable consequence of this scenario is an ex- 
periment in which a torsional oscillator containing solid "'He 
is set in motion, and finally is abruptly stopped and then let 
go again. If the sample is first cooled below the glass tran- 
sition while at rest, then the vessel contains a solid glass and 
the torsional oscillator will stop perfectly. Because all of the 
material inside the oscillator will stop in unison with the con- 
tainer. However, if the oscillator is stopped abruptly and let 
go at temperatures above the glass transition, then the vessel 
contains a liquidlike component and the torsional oscillator 
will wobble due to the motion of the residual liquidlike com- 
ponent. Here, the liquidlike component will still move even 
after the container is stopped. This internal Uquidlike motion 
will lead to a momentum transfer to the outer torsional oscil- 
lator once the container is let go. Hence, the oscillator will 
wobble. 



VIII. CONCLUSIONS 

In conclusion, we report on a simple alternate proposal of 
the origin of the decrease in the torsional oscillator period in 
solid heliumi 1 1- 12, 13. 14, 18. 19,20.21 argue that changes in the 
oscillation period are triggered by changes in damping and 
the effective stiffness in the susceptibility and are not due to 
a supersolid state formation. We point out that the torsional 
oscillator mechanical properties can be controlled by a num- 
ber of effects such as damping, stiffness coefficient and mo- 
ment of inertia changes. The supersolid interpretation, if such 
a state exists, would lead to changes in moment of inertia. 
This was the interpretation taken in Refs. [Tllfl2lll3lll4 On the 
other hand, low-temperature quenching of residual topologi- 
cal defects in the solidification into a glass can also account 
for the current observation of a drop in the torsional oscilla- 
tor period as the temperature is reduced. This effect not only 
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accounts for the decrease in the rotation period (Fig. [2]), but 
also predicts a peak in the dissipation, see Fig. [1] as indeed 
observed.— Furthermore, our ahernate glass scenario can ac- 
count for the observed dependence of the oscillator dynamics 
on external parameters. 

We have developed a basic understanding of how a simpli- 
fied response of a glass can account for the essential physics 
results of the torsional oscillator experiments. Specific heat 
data further show very strong support for a low-temperature 
glass phased and effectively rule out a simple supersolid tran- 
sition by thermodynamic measurements. Based on our anal- 
ysis, we can make the following specific predictions which 
conform with the data available to date: 

(1) At low temperatures all dynamic observables, e.g., Q^^ 
and P, will collapse on a single curve as a function of temper- 
ature, which is independent of the frequency of the torsional 
oscillator This collapse follows from Eqs. ( fTST i, ( fT9] l, and 
( |29] l. It reflects the freezing-out of liquidlike degrees of free- 
dom, where the response is controlled by the normal suscep- 
tibility xo ■ In general, the higher the resonant frequency ojq 
is, the higher the temperature T* at which this freeze out will 
occur: s(T*) = l/i^o, as can be seen from Eq.(|29]l, T* > Tq. 
At temperatures T < Tq, all of the data collected from oscil- 
lators with different resonant frequencies, in which all other 
parameters are held fixed, should collapse onto a single curve. 

(2) In all known glasses the exponent /? < 1 [Eq.(|20li1 and 
the glass temperature To > K, [Eq. (|29]l1, while in our anal- 
ysis we assumed for simplicity and to avoid the use of too 
many adjustable parameters that (3=1 and Tq = K. For 
more realistic glass models one would need to relax these con- 
straints. Indeed, Fig.([3]l vividly illustrates that better fits may 
be obtained with values of (3 < 1. Similarly, all of the avail- 
able data to date are consistent with values of < Tq < 100 
mK. The more accurate detailed predictions [albeit with more 
phenomenological glass parameters] follow from Eqs. (12111241 
|29l l. We leave this analysis for a separate discussion when 
more data are available. 

(3) The blockage of the annulus in the torsional oscilla- 
tor cell may drastically change the strain field and dynam- 
ics inside the cell. Such changes of the boundary conditions 
can inhibit the presence of various elastic modes relative to 
those which are present in a system with no blockage. For 
instance, azimuthal elastic modes which may be present in 
an "O-type" cylindrical geometry cannot appear in a blocked 
"C-type" cylindrical configuration. More generally, any mod- 
ification to the boundary conditions, e.g., oscillator velocity 
or amplitude, sample surface or roughness, would shift the 
freezing-out temperature and damping. For instance, large 
amplitude vibrations may thwart the quench of a liquidlike 
component into a glass. Such an avoidance of the glass transi- 
tion will, consequently, lessen the low temperature variations 
in the oscillator dynamics. All of these effects will necessarily 
lead to different experimental observables without a need for 
invoking a supersolid fraction which somehow depends very 
sensitively on all of these parameters. 

As a matter of principle, our glass picture explains the data 
even if a glass transition accompanies the transition into a 
supersolid statei^i^l What we show is that whether or not a 



tiny supersolid fraction is present- the changes in the dynam- 
ics can be naturally ascribed to a glass transition rather than 
to a Non Classical Rotational Inertia (NCRI) effect. Our re- 
sults demonstrate that a low-temperature glass can very natu- 
rally account for the current data. What our analysis suggests 
is whether or not a very small supersolid fraction is present, 
what triggers the change in period may (very naturally) be a 
hallmark of the transition into a glassy state. The changes in 
the torsional oscillator period need not be accounted for by 
only a non clasical rotational inertia effect. A glass transi- 
tion (whether it is classical or quantum is immaterial on this 
level of our analysis) may further naturally account for the ob- 
served subtle dependence of the effect on external parameters 
as well as the concurrency of the peak in the damping with the 
change in oscillator period. It is worth noting that a decrease 
in the rotational period along with a peak in the dissipation 
can be triggered by a host of many other related effects - none 
of which relies on a NCRI effect. As an example, in the ap- 
pendix we detail an alternate model in which the dissipation 
decreases by the solidification of liquidlike (topological) de- 
fects into a non-glassy low-temperature phase. Similarly, any 
other effect which leads to a decrease in the oscillator period 
with decreasing temperature and exhibits a peak in can 
account for the current data. 

Finally, an experiment in which a torsional oscillator hold- 
ing ''He is first cooled when it is at rest, then set in motion and 
abruptly stopped and let go, might allow further comparison 
with the glass origin that we propose here for the anomalous 
low temperature dynamics. For a glass transition scenario, 
we expect (1) that at temperatures above the glass transition, 
the liquidlike component will cause the torsional oscillator 
to wobble. While (2) the same experiment performed suffi- 
ciently below the glass transition temperature, the glass and 
container will stop in unison without wobbling. 
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APPENDIX A: AN ALTERNATE EFFECT - DISSIPATION 
LOWERING UPON FREEZING 

Here, we show how the seminal qualitative features of the 
current available dat a"''^i'^''^''^-'^'^'''^' can be explained by 
not only a glass transition [or even by its more restrictive lim- 
iting form, which is given by the simple activated dynamics 
of Eqns. (|20l [SOjl with (3=1 plotted in Figs. [T] and IS, but 
also by other, more general, "freezing" transitions. Glasses 
adhere to specific, nearly universal, response functions and 
relaxation time forms - e.g. that of Eq. (|20lii ^^-^^i^^ We now 
consider another "freezing" scenario of (topological) defects 
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in which, for illustrative purposes, only the dissipation varies 
with temperature. Upon cooling, the liquidlike regions in the 
solid "freeze" and lead to a lower average dissipation. 

In the dissipative response gdiss of Eq. (|8]l, the dissipation 
7(T) = jHe{T) + 'Yosc varies with temperature while /g/ / = 
losc+Ins is fixed. Thus, a measured value of the period can be 
used to compute ^{T) and then predict the amplitude damping 
rate Kdiss {T) and compare it to experiment. Here, 



Kdiss {T) 



lose + iHejT) 



leffP 



(Al) 



with lef f = lose + Ins HOW a coustaut quantity which is in- 
dependent of the temperature. Similarly, the Q factor (deter- 
mined by k) is 



Qdiss — 



(A2) 



If this alone is what occurs then at all temperatures, the dis- 
sipation rate k and Q-factor will be related to the period by 
Eqns. dAll IA2b parameterized by the two constants leff and 
a. On the right hand sides of Eqns. ( lAll IA2b . only the mea- 
sured oscillator period k varies with temperature, as indeed 
observed.'M^-'3.'4."^i'V°.^' 

The Q factor monotonically increases as the temperature 
decreases. This is evident from the final form of Eq. (IA2b or 
even from the trivial Eq. ( fT3] l. Because 7 decreases as the 
temperature is lowered, an ensuing increase of Q follows. 

We will now address the non-monotonic character of the 
dissipation. This effect is not present in the treatment thus 
far. The dissipation (probed by Q~^) was seeni^ to exhibit a 
pronounced maximum in nearly the temperature range where 
the oscillator period shows its most dramatic decrease. Now, 
near transition points (whether those of a defect freezing or 
of other origins), fluctuations and consequent dissipation are 
very often large. Here, we cannot think of only one damped 
oscillator for either the quenching of defects nor other tran- 
sition scenario. To make contact with the previous approach, 
we may divide the system into external torsional oscillators 
coupled to a multitude of small damped oscillators within the 
medium. If the medium is homogeneous then the system can 
be thought of as a single mode oscillator. Otherwise, a more 
careful analysis is required. In general, we will need to diag- 
onalize the N x N matrix [a — ijuj — ui'^I] with the elements 
lij, ctij and 7y being the parameters appearing in the equa- 
tions of motion for the angles \Bi\f^i of different elements of 
the solid. 

In the context of defect freezing, we might anticipate some 
defects to already be frozen (and lead to a lower value of 7 in 
those volume elements) and for other volume elements to have 
fewer frozen defects and to correspond to a higher value of 7. 
This distribution of damping coefficients leads to a solution 
for which is a superposition of many modes dispersed in 
frequency. The broader the frequency dispersion- the larger 
the incoherence- the smaller the value of the ((P') and that 



of the energy. As a result, the l/Q value which records the 
damping of the energy may be much higher near the transition. 

We next implement the stochastic character of the value of 
7 for each volume element. For illustrative purposes let us 
consider a distribution given by 



1 



2^a2 



:Cxp[-(7-7o)V(2a^)]. 



(A3) 



The width a-^ is generally a function of temperature. It is 
finite within the transition region and tends to zero outside 
it. Far from the transition, we have a single value of 7 [i.e., 
P(Ci) = ^(7 — 7o)] ™d our former analysis applies. Within 
the transition region 70 can be read off from the period by the 
use of Eq. (fTTT i. 

The angle is now a superposition of all modes. For a delta 
function torque Text{t) = A6{t), 

e{t)^Re{A d7P(7)6l^cxp[-iwJt-KTt]}, (A4) 
Jo 

g(j) the distribution of 7 values, which we anticipate to be 
approximated by Eq. ( IA3I ). and with uj^ and r'' given by 
Eqns. (EJinii. 

To get an intuitive feeling for the damping triggered by the 
dispersion, we can expand the argument in Eq. ( IA4I ) about 70 
to obtain the original solution (that assuming a single mode 
defined by a unique 70) multiplied by the approximate damp- 
ing factor 



C = 



Here, z = (7 — 70) and 



1 f"' 

/ dz ex.p[~iatz 

2TTa?, J -00 



zt 



21, 



eff 



(A5) 



(A6) 



Let us first examine an ideal situation where no damping ini- 
tially exists {k* = 0). Here, Eq. (IA5b is none other than the 
Fourier transform of a Gaussian which is another Gaussian 
given by a time dependence, C = cxp[— ia^cr^t^]. This fac- 
tor leads to an energy dissipation in time that does not appear 
for a single mode dispersion: the energy is damped by a factor 
of relative to that of the nondispersive system. Dispersion 
alone can lead to damping. Taking dissipation into account 
leads to a similar Gaussian factor of 



1 



C = cxp[--a:^i^(a^ + 



1 



4/2 



)] 



(A7) 



for the amplitude along with an additional small correction to 
the phase terms = uj^t uj^t + ^j^^^. This additional 
phase factor is negligible (over a period of the oscillator) if 
(4a/ -7^)3/2 > 2cr27. 

The above Gaussian approximation is just a heuristic moti- 
vation for visualizing how dispersion lowers the average en- 
ergy. The precise integral is that of Eq. ( IA4I ). If the width 
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(t(T) peaks at the transition temperature then the dissipation elude, we illustrated, once again, how the torsional oscillator 
is amplified near the transition temperature while, all along, data can be accounted for without invoking a supersolid tran- 
the oscillator period is monotonic in temperature. To con- sition. 
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